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Abstract 


We investigate how the criteria for two prototype instabilities in one-dimensional pattern- 
forming systems, namely for the Eckhaus instability and for the Benjamin- Feir instability, 
change as one goes from a continuous bifurcation, to a spatially periodic or spatially and/or 
time periodic state, to the corresponding weakly inverted, i.e. hysteretic, cases. We also 
give the generalization to two-dimensional patterns in systems with anisotropy as they 
arise for example for hydrodynamic instabilities in nematic liquid crystals. 

PACS classification: 47.20.Ky, 05.70.Ln, 47.25.Qv 
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1. Introduction 


The question of which wavelengths are linearly stable in a system showing spatially 
periodic states has been of interest since about 25 years dating back to the classical work 
of Eckhaus 1 . He investigated for a continuous bifurcation to a spatially periodic state 
how a given pattern with spatial variations in one direction (called one-dimensional pat- 
tern in the following) goes unstable upon applying a compressional or dilatational force. 
This instability, the Eckhaus instability, has since been studied theoretically for a num- 
ber of pattern-forming systems showing a forward bifurcation, including Rayleigh- Benard 
convection 2-4 and the Taylor instability 5 . Experimentally this one-dimensional instability 
has been studied first for the onset of electroconvection in nematic liquid crystals 6 and it 
has then been investigated in detail very carefully for Taylor vortex flow 7 . 

For the case of a continuous bifurcation to a spatially and time periodic state in 
one dimension the question of a modulational instability was addressed first for a purely 
dispersive system (i.e. no dissipation) by Benjamin and Feir 8 , when they investigated how 
water waves go unstable. Their analysis has been generalized by Newell 3 ' 9-11 to include 
both dissipation and dispersion. Later on in ref.12 this analysis was further extended by 
considering the wavenumber band possible above onset of the instability and not only the 
critical wavenumber as in refs. 3, 9-11. In addition, it has been shown 13 that the Benjamin- 
Feir-Newell instability is equivalent to the phase diffusion coefficient in the Kuramoto- 
Sivashinsky equation becoming negative. Up to the present there seems to exist no clear- 
cut experimental observation of the Benjamm-Feir-Newell instability in one-dimensional 
patterns. 

More recently there is an increasing interest in the spatial patterns formed near weakly 
inverted bifurcations (which show a small amount of hysteresis), especially when this bi- 
furcation is oscillatory in nature 14-20 . Motivated by these recent developments, which are 
partially triggered by experiments on thermal convection in binary fluid mixtures 21,22 , we 
investigate here how the criteria for the Eckhaus instability and the Benjamin-Feir-Newell 
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instability change as one goes from a forward to a weakly inverted bifurcation. 


2. The Eckhaus instability for a weakly inverted bifurcation 
to a spatially periodic state 

The envelope equation for a weakly inverted bifurcation to a spatially periodic state 
reads 

A — eA + yA xx +(3\A\ 2 A - 8\A\ A A ( 2 . 1 ) 


where e, 7 > (3 and 6 are real and where f3 and 6 are assumed to be positive to guarantee 
that the bifurcation is inverted and saturates to quintic order. The diffusion coefficient 7 
will also be assumed to be positive. 

The time-independent version of eq.( 2 . 1 ) allows finite amplitude, plane wave solutions 
of the form 

A at = A 0 e ikx ( 2 . 2 ) 


where 




(2.3) 


= W) 1 1 ± 

In eq.(2.3) the upper (+) sign in the bracket refers to the physical branch of solutions, 
which exists for a band of k values 


7 k 2 < 



(2.4) 


whereas the lower (-) sign is associated with the unphysical branch of solutions. 

To test for linear stability of the stationary solutions A st (eq.( 2 . 2 )) against spatial 
modulations, we proceed along the lines of ref .2 and write 


A — A at -f Ai(x , t ) 

We obtain upon linearization in A\ 


(2.5) 


i, = eA\ +jA lxx - A 1 \A 0 \ 2 (-2(3 + M\A 0 \ 2 )-A 2 0 A* 1 e 2ik *(-{3 + 28\A 0 \ 2 ) ( 2 . 6 ) 
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This equation allows for solutions of the form 


Ai(z,t) = ^u(<)e i(fc+M) * + A i2 (t)e^ k ~ M ^ 


(2.7) 


Setting An, A 12 oc e M gives 


A 1)2 = -( 7 (k 2 + M 2 )-e + \A 0 \ 2 (-2(3 + 35|A„| 2 )) ± + |A 0 | 4 (-/? + 25|^ 0 | 2 ) 2 

( 2 . 8 ) 


A 1]2 (M,fc) has a relative maximum for 


M 2 = k 2 


|Aoj 4 (-/3 + 2j»l^o| 2 ) 2 
4 7 2 fc 2 


(2.9) 


and 

A = e - \Ao\ 2 (-2(3 + 35|A 0 | 2 ) + ^(-/? + 25|4„| 2 ) 2 (2.10) 

In deriving eq. (2.10), the positive sign was chosen in eq. (2.8) since it corresponds to the 
larger growth rate. It is easily checked from eqs. (2.1), (2.9) and (2.10) that for 5 = 0 and 
j3 negative the classical result for a forward bifurcation 1 ’ 2 (which reads 


Jfc 2 < 


37 


( 2 . 11 ) 


in the present notation) is obtained. Inserting eq.(2.3) into eqs.(2.9), (2.10) we find that 
the finite amplitude solutions (eqs. (2. 2, 2.3) with positive sign) of eq.(2.1) are Eckhaus 
stable, i.e. stable against the growth of sidebands, provided 


k 2 < 


f 3/3 2 
2 7 32 7 5 



( 2 . 12 ) 


Inequality (2.12) represents the stability domain for stable finite amplitude plane wave 
solutions of eq.(2.1) that are stable against the Eckhaus sideband instability. Larger values 
of the wavevector will lead to an instability. Comparison of inequalities (2.12) and (2.4) 
shows that the Eckhaus stable range is smaller than the existence range of nonlinear plane 
wave solutions of the type (2.2), (2.3). In analogy with the experimental observations 6 ’ 7 
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and with theoretical investigations 23 for a forward bifurcation one can expect the gain or 
loss of roll or vortex pairs to occur for the Eckhaus instability. This expectation is based on 
the fact that one has the same mechanism (a sideband instability) for both cases, namely 
for the spatially periodic states formed for a forward and for a weakly inverted bifurcation. 

3. The analog of the Benjamin-Feir-Newell instability 

For a weakly inverted bifurcation to a spatially and time periodic state one has for 
spatial variations in one dimension the envelope equation 13-2 ®’ 2 ^’ 23 

A = eA + iA xx - f3\A\ 2 A - S\A\ 4 A (3.1) 

where 7, /?, and 6 are complex and thus of the form a = a r + ieti and where fi r < 0 
and 8 r > 0 to guarantee that the bifurcation is inverted and saturates to quintic order. 
Furthermore we have discarded nonlinear gradient terms in writing down eq.(3.1). 

Eq.(3.1) admits spatially homogeneous solutions of finite amplitude 


Ah — Aoe 1 


where 


w — — Mo| 2 (/?t + £t|-<4o| 2 ) (3.4) 

and where the (+) sign in the bracket of eq.(3.3) corresponds to the physical branch, 
whereas the (-) sign is associated with the unphysical branch. In the spirit of refs. 9-11 
we investigate the stability of the spatially homogeneous solutions (eqs.(3.2)-(3.4)) with 
respect to space-dependent perturbations, that is we look for solutions of eq.(3.1) of the 


A(x,t) = A h (t) + B(x,t) 
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Inserting eqs.(3.2) - (3.5) into (3.1) we obtain upon linearization in B 

B = eB + 7 B xx - B\A Q \ 2 {2j3 + 36|A 0 | 2 ) - A 2 B*e 2iut {(3 + 28\A 0 \ 2 ) (3.6) 

Using the ansatz 

B = Bi(i)e ikx+iu>lt + B 2 (t)e~ ikx+iUit (3.7) 

and 2u; = u>\ + u) 2 one gets a closed system of equations for Bi(t ) and B%(t). Setting 
i?i(<) oc B ioe At and B 2 (t) cx B 2 oe Xt , we arrive at the following solvability condition for A 

Al ( 2 = — {jrk 2 ~ t + 2/? r |j4o | 2 + 38 r \Aq | 4 ) 

± VlMWr + 2^r|-do | 2 ) 2 + Wi + 2^t|j4o| 2 ) 2 ) - ( 7 i* 2 +o; + |Ao| 2 (2 ft + 3^|A 0 | 2 )) 2 

(3.8) 

Upon replacing t and u> in eq. (3.8) with the values given by eqs. (3.3) and (3.4), re- 
spectively, analysis of eq.(3.8) shows that the spatially homogeneous solution (3.2) - (3.4) 
is linearly unstable against spatially inhomogeneous perturbations of the type considered 
here if 

7t {.Pi + 2<$i|>lo| 2 ) + 7 r(/?r + 28 t \Aq\ 2 ) < 0 (3*9) 

Eq.(3.9) represents the analog of the Benjamin-Feir-Newell criterion for a weakly inverted 
bifurcation. For 8 r = Si = 0 the case of the forward bifurcation 3 t9 " 11 is easily recovered 
from eq.(3.9) 

Jifii + Irftv < 0 (3.10) 

As for the case of the forward bifurcation 9 " 11 , the terms in eq.(3.9) represent a competition 
between the ’cooperative’ (dissipative and diffusive) tendencies to generate an ordered 
structure at q = go (the wavevector of the underlying structure such as for example rolls 
or vortices) and the dispersive forces (dispersion and nonlinear refractive index), which try 
to focus initial nonlinear wavetrains into pulses. 

To investigate the modulational instability of plane wave solutions with a finite 
wavenumber we proceed along the lines of ref.12. Eq.(3.1) has finite amplitude plane 
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wave solutions of finite amplitude of the form 


A p 


Aoe i (kx+urt) 


(3.11) 


where 


and 



Pt_ 
2 6 r 




= -H k2 ~ \ A o\ 2 {Pi + ^iMol 2 ) 


(3.12) 


(3.13) 


and where the (+) sign in the bracket of eq.(3.12) corresponds to the physical branch, 
whereas the (-) sign is associated with the unphysical branch. To investigate the linear 
stability of the finite amplitude plane wave solutions (eqs.(3.11)-(3.13)) with respect to 
space-dependent perturbations is we look for solutions of eq.(3.1) of the form 


A(x,t) = A p (x,t) + B(x,t) (3.14) 

We find that, linearizing in B, structurally the same equation for B is obtained (eq.(3.6)) 
as for the case of the spatially homogeneous solutions, except that now \A 0 \ 2 and lo are 
given by eqs.(3.12),(3.13) instead of eqs.(3.3),(3.4). 

It is possible to find a solution of the resulting equation for B of the form 

B(x,t) = a(t)e i(klX+u>lt) + b(t)e i(k **+“»*) (3.15) 

provided k\ + k 2 = 2 k and uq -(- u >2 = 2u>. Introducing the abbreviations 


ot r — c + 2/? r |j4o| 2 T 3^ r |^4 0 | 4 

(3.16) 

- 2/3 i |^o| 2 + 3^ i |^o| 4 

(3.17) 

< r =0 r + 28 r \ A »\ 2 

(3.18) 

<<=& + 26i\A 0 \ 2 

(3.19) 
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we find for a(t) and b*(t) the equations 


a = 

— ((7r + *7«) fc l + ttOl + («r + icti))a - (Cr + iC*)A lb* 

(3.20) 

b* = 

((TV ili) k 2 *^2 "f* (c*r (C* - *C«)^0 ® 

(3.21) 

Assuming the disturbances to be proportional to e At , we find that A has a negative real 

part provided that 

— — (Si + £ 2 ) — OL r + Ncos(~9) < 0 
z z 

(3.22) 

where 




1 1 \ * 
■Si) 2 - 4 (r, + r 2 ) 2 + A 2 ) 2 + -(r, + r 2 ) 2 (s, - s 2 ) 2 J 

(3.23) 

and 

w =1 5 < r > + r 2 )(Si -s 2 ) 

|(s 1 -s 2 )2-|(r n +r 2 )2 + A2 

(3.24) 

with 

r j — cti + jtkj + a 

(3.25) 


S . j = Irk) 

(3.26) 


> 

II 

lu 

0 

_ 

(3.27) 


where j = 1,2. Comparing eqs.(3.22)-(3.27) with eqs.(3.14) - (3.18) of ref.12, we see that 
there is a great similarity in structure. Further analysis of inequality (3.22) turns out 
to be rather involved for the case of a weakly inverted bifurcation and is left to future 
investigations. In closing this section we note, that eqs.(3.22) - (3.27) contain as a special 
case inequality (3.9), as one sees after some algebraic manipulations. 

4. Generalization to anisotropic systems 

So far we have looked at envelope equations containing only spatial variations in one 
spatial dimension. Now we generalize this to the case of spatial variations in two dimensions 
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for systems which have an intrinsic (such as anisotropic liquids: nematic liquid crystals 
etc.) or extrinsic (e.g. due to an external force) anisotropy as it occurs for example for 
hydrodynamic instabilities in magnetic liquids in an external magnetic field. 

In this case the envelope equation (2.1) for a weakly inverted bifurcation to a spatially 
periodic state is replaced by 

A = eA + 711 A xx + 712^ + 722 Ayy + ( 3 \A\ 2 A - < 5 |A| 4 A ( 4 . 1 ) 


with 7n, 712, and 722 real. 

An envelope equation for a forward bifurcation to a spatially periodic state has been 
derived for electroconvection in nematic liquid crystals 26 . As for the case of the forward 
bifurcation 26 one can diagonalize the terms containing second derivatives of the envelope 
in eq.( 4 . 1 ) by rotating the coordinate system and by rescaling the coordinates. Then 
the derivative terms take the form of a two-dimensional Laplacian A2 and the resulting 
envelope equation for a weakly inverted bifurcation to a spatially periodic state reads 

A = eA + 7 A 2 A + fi\A\ 2 A - S\A\ 4 A (4.2) 


It is easy to show that the most unstable perturbations correspond to perturbations along 
the plane wave and therefore all the steps we have outlined in section 2 for the Eckhaus 
instability in one dimension apply. We therefore obtain for the band of Eckhaus-stable 
wavevectors 


k l+ k l< o - + 


3 /? 2 


1 + 4/1 + 


32 eb 
9 ( 3 2 


( 4 . 3 ) 


27 327£ 

where £ and 77 in inequality ( 4 . 3 ) denote the coordinates in which the spatial derivative 
terms in eq.( 4 . 1 ) assume the structure of the two-dimensional Laplacian. 

For a weakly inverted bifurcation to a spatially and time periodic state in an 
anisotropic system, eq.(3.1) is replaced by 


A - eA + ~i\\A xx + 712 A X y + 722 Ayy - P\A\ 2 A - 8 \A\ a A ( 4 . 4 ) 
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where 7n, 712, and 722 are complex and of the form 'fkl — 7 + iYki- For the envelope 
equation (4.4) we can investigate the analog of the Benjamin-Feir-Newell instability along 
the same lines as for the one-dimensional case in section 3. We find that the analysis goes 
through and that inequality (3.9) is replaced by the requirement that the (2x2) matrix Efcj 

Sfci = 7 li(0i + 2^i|^4o| 2 ) + Ikiiflr + 2S r \A 0 \ 2 ) (4-5) 

must be positive to guarantee linear stability of finite amplitude plane wave solutions 
against space-dependent perturbations. 

This analysis might turn out to be useful in the near future, since an inverted bifur- 
cation to traveling waves has been found recently experimentally for the onset of electro- 
convection in nematic liquid crystals 27 . In addition it has been predicted a few years ago 
theoretically 28 , that the onset of thermal convection in nematics can also arise as a weakly 
inverted oscillatory instability. We note that this calculation has been done for homeotropic 
alignment of the director. The preferred direction in the plane is provided by the external 
destabilizing magnetic field, which is applied in a direction parallel to the nematic layer. 
Thus the planar horizontal symmetry is not broken by the orientation of the director due 
to surface treatment as in the planar configuration studied in electroconvection, but broken 
externally by the direction of the external in-plane magnetic field. 

For the case of a forward oscillatory bifurcation, the condition that the matrix in 
eq.(4.5) must be positive for stability is replaced by the requirement that the (2x2) matrix 

Sfci 

EfcJ — IklPi + IklPr (4.6) 

must be positive for stability, which is a direct application of a previous general theoretical 
result by Newell 10 to a real physical system, namely electrohydrodynamic instabilities in 
nematic liquid crystals. Envelope equations for a forward bifurcation to traveling waves 
have been discussed very recently 29 in an analysis of experimental results on electroconvec- 
tion in nematics 30 , where a transition from traveling to standing waves could be achieved 
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by modulation. 


5. Summary and perspective 

In this note we have discussed the analogues of the Eckhaus and the Benjamin-Feir- 
Newell instability for a weakly inverted bifurcation. In addition, we have generalized these 
criteria to the case of anisotropic liquids such as nematic liquid crystals. 

From an experimental point of view the next step could be the investigation of the 
Benjamin-Feir-Newell instability for traveling waves in thermal convection in binary fluid 
mixtures, which are well documented for the upper (stable) branch of the hysteresis loop. 
This will allow a test of the criterion given here, provided one choses a value for the 
separation ratio such that the size of the hysteresis loop is sufficiently small and thus an 
envelope equation is applicable. To test the criterion given for the Eckhaus instability for 
a weakly inverted bifurcation to a spatially periodic pattern one could also study binary 
convective mixtures. For such mixtures in a porous medium a regime of separation ratios 
has been found over which the bifurcation is expected to be inverted 31 ’ 14 . For the case 
of anisotropic fluids such as nematic liquid crystals, it will be most interesting to see 
under what circumstances the bifurcation to the stationary branch is actually forward or 
inverted. In the latter case, which has been observed for the onset of electroconvection 32 , 
one could test the criterion presented here for the Eckhaus instability near a weakly inverted 
bifurcation to a stationary pattern. 


Note added-. After submitting this paper for publication, we became aware of ref.33. 
This paper contains in eq.(1.5) implicitly our criterion (eq.(3.9)) for the analogue of the 
Benjamin-Feir-Newell instability for a weakly inverted bifurcation. Due to the scaled form 
of eq.(l.l) in ref.33, however, it is not possible to go through the tricritical point from 
the weakly inverted case to a forward bifurcation. In our description this can be done 
straightforwardly (compare section 3). 
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